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Abstract. We study a Dirac operator subject to Atiayh-Patodi-Singer like boundary con- 
ditions on the solid torus and show that the corresponding boundary value problem is elliptic, 
in the sense that the Dirac operator has a compact parametrix. 



1. Introduction 

The celebrated Atiyah-Patodi-Singer boundary condition [2] for Dirac operators on closed 
manifolds with boundary was introduced as a key ingredient in the generalization of the 
Atiyah-Singer index theorem. It is a non-local boundary condition which makes the Dirac 
operator Fredholm. An advantage of the APS condition is that further detailed analysis 
can be carried out where local conditions such as Dirichlet and Neumann may not be well 
behaved. 

That theory works, however, under the assumption that both the manifold and the op- 
erator have a cylindrical structure near the boundary. A semi-infinite cylinder can then be 
smoothly attached to the manifold and the Dirac operator can be naturally extended over 
to the cylinder. The APS boundary condition can be then described in the following geo- 
metrical terms: a sufficiently regular section is in the domain of the operator if it extends to 
a square integrable solution on the cylinder. The issue though is that many concrete nat- 
ural operators do not have a cylindrical structure near the boundary. Even a simple d-bar 
operator d/dz on a disk in the complex plane does not have this structure. The solid torus 
example studied here is not cylindrical near boundary either. 

In this paper we consider a Dirac operator on the solid torus considered geometrically as 
the product of the unit disk and the unit circle. We construct another non-local boundary 
condition similar in spirit to the APS boundary condition. It was inspired by the non- 
local boundary conditions discussed in [12] and [6] to get around the necessity of having a 
cylindrical structure near the boundary. The boundary condition we propose in this paper 
has the same geometrical interpretation as the APS condition. Namely, we consider the solid 
torus as a subset of the bigger noncompact space of the plane cross the unit circle. We define 
the domain of the Dirac operator in full analogy with APS as consisting of those sufficiently 
regular sections which extend to square integrable solutions on the complement of the solid 
torus. 

The motivation for studying this particular example comes from our larger project of devel- 
oping a concept of a noncommutative manifold with boundary and noncommutative elliptic 
boundary conditions. This is done by studying examples, starting with two-dimensional do- 
mains and continuing with more complex cases. In particular our efforts were concentrated 
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on studying quantum analogs of Dirac operators subject to APS like boundary conditions, 
see [5], [8], [9], and [10]. The solid torus studied in this paper is possibly the simplest 
three dimensional example, yet significantly more difficult then the two-dimensional exam- 
ples studied in our previous papers. While the standard APS theory does not apply, the 
example however seems to have an attractive noncommutative version. In the spirit of the 
series, this paper will be followed up by another study that discusses the same type of theory 
and calculations for the quantum analog to the solid torus. 

We show in this paper that the Dirac operator on the solid torus subject to our non- 
local boundary condition is self-adjoint and has no kernel. Using a partial Fourier transform 
we obtain an explicit formula for its inverse. We then show that the inverse is a compact 
operator, which is the main feature of elliptic boundary conditions. Moreover we prove that 
the inverse is a p-th Schatten class operator for p > 3. This is obtained by direct analysis of 
the formula for the inverse using subtle estimates involving modified Bessel functions. 

The paper is organized as follows. In section 2 the Hilbert space and the Dirac operator D 
are defined, and the boundary condition is stated. The section also contains the computation 
of the kernel of the Dirac operator, and the computation of its inverse Q. The proof of the 
main theorem, the compactness of Q, and the Schatten class computation is contained in 
section 3. In the appendix we collect numerous facts about the modified Bessel functions, 
some are classical and some are more recent. 



2. The Dirac Operator 

We begin with the necessary notation. Let D = {z G C : \z\ < 1} be the unit disk and let 
S* 1 = \e %e G C : < 9 < 27r} be the unit circle. Also let T 2 be the two dimensional torus and 
let ST 2 be the solid torus: ST 2 =3xS 1 C Cx S 1 . The boundary of ST 2 is T 2 . The operators 
that we are studying will be acting in the Hilbert space U = L 2 (ST 2 , C 2 ) = L 2 (ST 2 ) <g> C 2 , 
i.e. the space of square-integrable complex vector-valued functions on the solid torus. The 
inner product of F, G G V. will be denoted as (F, G) . 

We proceed to the definitions of the main object we study in this paper. We consider the 
following formally self-adjoint Dirac operator D defined on % by 

/ ld_ r)d_ \ 

V dz i 80 J 

Notice that D can act on either L 2 (l x 5 1 ) ® C 2 or L 2 (C x S 1 ) <g> C 2 . 
The domain of D is defined to be: 

dom(D) = {F e H\ST 2 ) ® C 2 : 3F ext G H] oc ((C x S 1 ) \ ST 2 ) ® C 2 } (2.2) 

such that (1), (2), and (3) hold: 

(1) F £Xt \ T 2 =F\j2, 

(2) DF ext = 0, 

(3) F ext G L 2 ((C x S 1 ) \ ST 2 ) <g) C 2 . 

Here H l is the first Sobolev space. 

The first task is to study the kernel of D. For a function F G L 2 (ST 2 ) <g> C 2 , by using the 
polar decomposition z = re llp , we have the Fourier decomposition as follows: 
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The norm of F can then be expressed as 



p _ \ r I fm,n( r ) \ ^inip+ime ^ 3) 



\\F\\ 2 =(F,F) = f (\fm,n\ 2 + \gm,n\ 2 )rdT. 

m,neZ 

We wish to solve the equation DF = without any kind of conditions imposed. This is 
done in the following proposition. 

Proposition 2.1. Let D be the operator defined by (2.1) and acting in the Hilbert space 
L 2 (ST 2 \ ({0} x S 1 )) <g> C 2 . Then the kernel of D consists of those F e L 2 (ST 2 \ ({0} x 
S 1 )) <g> C 2 for which the coefficients of (2.3) satisfy the following relations: for m ^ and 
any n 

m 

fm,n+i(r) = | — | (-A mi „7 n+ i(|m|r) + B m , n K n+1 (\m\r)) (2.4) 

and 

), (2-5) 

while ifm = and any n then f , n (r) = A 0)Tl r~ n and g , n (r) = B 0)Tl r n . Here A m ^,B m ^ n are 
constants and I n , K n are the modified Bessel functions of the first and second kind respec- 
tively. 

Proof. To solve the equation DF = 0, it is enough to solve the equation 

Ef m e%L<> (J^ ~ V £)\( fm,n(r) \ in<p+ime _ n 



by the Fourier decomposition (2.3). There are two cases to consider. First for m ^ and 
any n, letting t = \m\r, f m , n (t) = f m , n (t/\m\), and g m , n {t) = g m , n (t/\m\) we arrive at the 
following system of differential equations: 

~ n ~\~ 1 ~ m 

ClW + —r-fm,n + l(t) + T-r~g m ,n(t) = 
6 

Tl 771 ~ 

g'm,nif) - -jgmAt) + j^r/m,n+l(*) = °" 

Substituting the second equation into the first yields: 



~C, n {t) + ^ - (l + £ ) g m , n (t) = 0, 



which is equation (A 3). This implies that g m , n (t) is a linear combination of the two modified 
Bessel functions, i.e. g m>n (t) = A m>n I n (t) + B m n K n (t) . Then, using the second equation from 
the above system and using the relations (A. 4), we get the desired result. 
If m = then the above system of differential equations reduces to: 



4 



SLAWOMIR KLIMEK AND MATT MCBRIDE 



/o» + ~ r kn(r) = 
9o, n ( r ) ~ -9oA r ) = 0. 

This system is uncoupled and the solutions are easily seen to be /o,n( r ) — Ao )H r~ n and 
9o,n( r ) — B n r n . Thus this completes the proof. □ 

In the following proposition we rephrase the domain condition of D in terms of the Fourier 
coefficients to explicitly write down the boundary condition. 

Proposition 2.2. Suppose that F is in the domain of D and has Fourier decomposition 
given by (2.3). Then 

\m\K n+1 (\m\)g„ hn (l) - mK n (\m\)f mj „ +1 (l) = (2.6) 
ifm^O and fo, n (l) = for n < 0, and go, n (l) = 0forn>0. 



Proof. We want to find a function F ext for r > 1 such that DF ext = and F ext \ r= i = F\ 
First solving the system DF ext = for F ext yields 



r=l- 



fext t r \ 
J m,n\ J 

\ 

m.ngS 



E 



7<ext \ ( J m,n 

^ext 



~eXt { tr.\ 

ym,n\ J 



with fmn+i{ r ) an d J™n( r ) given by (2.4) and (2.5) for m^O. Additionally /o^*(r) = Aq^ 71 
and g^ni^) = B 0tTl r n . First consider m ^ 0. In order for the solutions to agree on the 
boundary of the disk, the coefficients A m , n and B m<n must solve the following system of 
equations 

I n (\m\) K n (\m\) \ ( A m , n \_f g m , n (l) 

-R^l(l^l) W\ K n + l(H) ) \ B m , n ) - \ f m>n+1 (l) 

The solution is: 

A m ,n = \m\K n+1 (\m\)g mjn (l) - mK n (\m\)f m:Tl+1 (l) 

and 

B m , n = \m\I n+1 (\m\)g„ hn (l) + ml n (\m\)f rn . n+ i(l) 
If m = we get /o, n (l) = A , n and g 0t n(l) = B 0)fl . 

Our boundary condition requires that F ext is square integrable on the complement of ST 2 
in C x S 1 . Because of the asymptotic properties of the modified Bessel functions, see (A. 8), 
this forces A mtTl = for m ^ 0. If m — the integrability of powers of r force A 0>n = for 
n < 0, and B 0:Tl = for n > 0. The statement follows from the above formulas for A m>n , 
A 0i „, and B 0:Tl . □ 

Combining the above two propositions we obtain the following corollary. 

Corollary 2.3. Let D be the operator defined by (2.1) subject to the boundary conditions 
(2.2). Then its kernel is trivial. 
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Proof. Let F e L 2 (ST 2 \ ({0} x S 1 )) <g> C 2 be a solution of DF = 0, as in Proposition 2.1. 
Its extension to a solution on C x S 1 \ ({0} x S* 1 ) is clearly given by the same formula. Since 
the powers of r are either not square integrable at zero or at infinity, it is clear that A n = 
and B 0tTl = 0. 

If m 7^ then in order for F to be regular at zero we must have B m n = by the asymptotic 
expansion of K n (t) near zero, see (A. 7). Then Proposition 2.2 implies that A m>n = 0, hence 
the kernel of D is trivial. □ 

It turns out that our boundary condition is self-adjoint as demonstrated in the next propo- 
sition. 

Proposition 2.4. The operator D defined by (2.1) subject to the boundary conditions (2.2) 
is self- adjoint. 

Proof. It is clear that D is formally self-adjoint. From the standard elliptic theory [4] the 
domain of D and its adjoint consists of (vector- valued) functions of the first Sobolev class. 
Thus the only thing that we need to check are the boundary conditions of the adjoint. To 
this end we inspect the boundary integral in Green's formula. Let F, G be H l functions on 
the solid torus. Using Proposition 2.2 from [8] and the Fourier decompositions: 

f \ \ r / /m,n(^) 



^inip+irad 



e 



we obtain: 

(DG, F) - (G, DF) = 2 «p, dg/dz) - (q, df/dz) - (dq/dz, f) + (dp/dz, g)) 

= 2 ^ (P™,n+l(^)gni,n( l ) ~ ?m,n(l)/m,n+l(l)) • 

Now suppose that F is in the domain of D, so it satisfies the conditions of Proposition 2.2. 
For G to be in the domain of the adjoint of D we need the above expression to be equal to 
zero. This gives: 

E( m K n (\m\) t— r-A „ 

P0,n+1 (l) 50 ,n(l)-^*,n(l)/0,n+l(l) =0. 

n<0 n>0 

For m^0 the above equation will equal zero for arbitrary F only if |m|ii' ri+ i(|m|)g miri (l) — 
mK n (\m\)p mtn+ i(l) = 0. Additionally we must have Po,n(l) = for n < 0, and go,n(l) = 
for n > 0. All together those requirements are exactly the same as the original boundary 
condition. Hence D and D* have the same domain and the proof is complete. □ 

The next goal is to construct the inverse of D. This is done by explicit solving of a non- 
homogeneous system of differential equations for the Fourier components and adjusting the 
integration constants to get the regularity at r = and so that the boundary condition is 
satisfied. 
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Proposition 2.5. Let D be the operator defined by (2.1) subject to the boundary conditions 
(2.2). Then the operator Q given by the formula (2.9) below is the inverse to D, in other 
words QD = DQ = I. 

Proof. To compute the inverse of D we solve the equation DF = G, which we will reduce 
to solving a non-homogeneous second order ordinary differential equation. We will use 
the Fourier decomposition (2.3). We first consider the case m ^ 0. Letting t = \m\r, 
fm,n(t) = f m , n (t/\m\), and similarly for other functions, DF = G becomes the system of 
differential equations: 

t \m\ \m\ 

f, ^ n + 1 f (f] m ~ (f) _ Qm,n(t) 

- Jm,n+lK T ) 7 Jm,n+l{t) ~ J— T9m,n{ T ) ~ i i • 

6 

Next by substituting the first equation into the second we arrive at the following second 
order differential equation: 



1 / / s ( . TV 



,2 



~€,n(t) + f9m,n(t) ~ (I + -p) 9m,n(t) = 

~i i \ ~ (2-8) 

J_~ ,.s , Pm,n+l{t) U+1 p m ,n+l(t) 

m \m\ t \m\ 

where the right hand side of the above equation was denoted by h m>n (t). Notice that equation 
(2.8) is the non-homogeneous version of the Bessel differential equation (A 3). General theory 
of ordinary differential equations tells us that the general solution to equation (2.8) is 

~9m,n(t) = ^(tWt) + cj*^ \t)K n (t) 

where cf 1 '^ (t) , c^ 1 ^ (t) solve the following system 

i n (t) K n ( t )\ ( c { r n \t)\ _( 

I'n(t) Kit)) { C ^ n \t)J ~ \h m , n (t) 

The solution of this system is 
cj m ' n) (t) = A m , n + / sK n (s)h mtn (s)ds and c { ™' n) {t) = B m)H - sl n (s)h m , n (s)ds, 

J\m\ Jo 

where A mjTl and B m<n are constants. The boundary condition (2.6) and the regularity at 
t — imply that we must have cf 1 '^ (t) equal to zero on the boundary, in other words where 
t = \m\. The boundary condition and regularity also imply that c^™'™^) goes to zero as 
t — > 0. After performing the integration by parts we obtain that A mjTl = K n (\m\)p min+ i(\m\) 
and B mtn = 0, and 



\m\ 

r-t 



c ^ n \t) = -—f S l n ( s )q ( s )d s + -L [ sl n+ i(s)p mtn+1 (s)ds. 
m Jo \m\ Jo 



1 


/' 


\m\ , 


l\m 


I 


f 


\m\ , 


JO 
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Next we use the first equation of (2.7) to solve for the general f m ,n+i(t) term to get: 
m f (+\ — Pm,n+i(t) n_ 

T—Jm,n+l(t) - r-| h -gmA 1 ) ~ 9m,n\P)- 

A straightforward calculation using the relations between the derivatives and indices of the 
modified Bessel functions (A. 4) and the formula for the Wronskian of the modified Bessel 
functions (A. 6) yields: 

fm,n + l(t) = ^- (-C { r' n \t)I n+1 (t)+ct' n \t)K n+1 (t)) . 

| TO I V / 

Next consider the case when to = 0. The system of differential equations reduces to the 
following 

9oA r ) ~ ~9oA r ) = Po,n+i(r) 

fo,n+i( r ) + ^y-hn+iir) = -q , n ( r ) 

This system is an uncoupled system and can be solved using an integration factor in each 
equation. Therefore the formula for the parametrix to D is 

QG := V ( f m ' n ^\ \ e inf+ime + V- f hn(r) \ j nlp 

where for to ^ 0: 

fm,n+l(f) — 

= \m\I n+1 (\m\r) J K n {\m\p)q m , n (p)pdp + ml n+1 (\m\r) J K n+l (\m\p)p m ^ n+1 (p)pdp 

- \m\K n+1 (\m\r) / I n (\m\p)q m>n (p)pdp + mK n+1 (\m\r) / I n+1 (\m\p)p m , n+1 (p)pdp 
Jo Jo 

9m,n{ r ) = 

= -ml n (\m\r) / K n (\m\p)q m , n {p)pdp - \m\I n (\m\r) j K n+1 (\m\p)p m<n+1 (p)pdp 

J r J r 

pr pr 

-mK n (\m\r) / I n (\m\p)q mtn (p)pdp+\m\K n (\m\r) / I n+ i(\m\p)p mjn+1 (p)pdp 
Jo Jo 



and 

fo,n+i(r) 

and 

9oA r ) = 



r p n 

jqo, n (p)p d P n > 



1 



P 

qo,n{p)P d P n<0, 



/ -^Po,n+i(p)pdp n>0 

J r P 

-^Po,n+i(p)pdp n<0. 
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It is is now a routine exercise to verify that DQ = QD = I. Thus this completes the 
proof. □ 

3. The Parametrix 

Now that the parametrix has been constructed the next goal is to show that it is a compact 
operator. This is the main result of the paper. 

Theorem 3.1. The Dirac operator D , defined by (2.1) and subject to the boundary conditions 
(2.2) has a bounded inverse. Moreover that inverse is a compact operator. 

Proof. Consider the following integral operators in L 2 ([0, l],rdr) for i, j = 0, 1 and m^fl: 
R ( ip n) f{r) := \ m \ J I n+i(\m\r)K n+j (\m\p)f(p)pdp, 

Sip n) f(r) := \m\ f K n+i (\m\r)I n+j (\m\p)f(p)pdp, 



and for n > 0: 



i(0,n) 



Jr P 



r 



Tt n) f(r):= / ^Mpdp. 
Jo ' 

Then we can rewrite formula (2.9) for Q in the following way: 

f - D( m ' n ) \ m \ D( m ' n ) o( m ' n ) I \ m \ q(m,n) 

Jm,n+1 — -n-10 Qm,n ~T -tin Pm,n+1 °10 Qm,n ~T On P m,n+1 

•" , ; "' (3.1) 

\m\ jj(rn,n) T->(m,n) \m\ ^( m ,n) Q(m,n) 

9m,n = -^00 Qm.n ~ -^01 Pm,n+1 ^00 Qm.n + ^11 Pm,n+1 

m m 



and 



and 



/o,n+l — 



-T 2 (0 ' n) g 0i „ n > 
T^~ n ~ l) q,, n n<0, 



90,n 



-Ti°' n) p 0tn+1 n>0 
T 2 (0 ' _n_1) po,n+i n < 0. 



We will show that all ten integral operators above are Hilbert-Schmidt by estimating the 
Hilbert-Schmidt norms. It turns out that the HS norm of each integral operator goes to 
zero as \m\ + |n| goes to infinity. This implies that Q is a compact operator as the norm 
limit of compact operators, since it is (up to a shift in the n index) a direct sum of compact 
operators with decreasing norms. 

To show that T^ '" - -' and 

T (0,n) 

are Hilbert-Schmidt we simply compute: 

For the other operators we have: 
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1 1 2 = I m 1 2 / / J n+i ( I m \ r ) K l+j ( I m I P) rpdpdr, 

JO Jr 



and 



11^112 = H 2 / f Kl +l (\m\r)ll +3 (\m\p)rpdpdr, 
Jo Jo 

Clearly we have R^' n ^ = i?Q™' n+1 " ) and S^' n ^ = Sq™'™^ and additionally, using the inequal- 
ity (A. 9), we can conclude that: 

||p(m,n)i|2 ^ ||p("i>n)||2 ||p( m > n )||2 
H-n-lO \\2 — 1 1 -"DO 112 — 1 1 -"-01 1 1 2 

I. q(m,n) II 2 ^ || q(m,n) m2 ^ ||c( m ' n )||2 
ll°01 112 — ll D 00 1 1 2 — 11° 10 1 1 2 - 

Consequently, we only have to estimate the Hilbert-Schmidt norm fori2<?' n) and S[^ n) . We 
have: 

1 r\ m \ r\ m \ 1 r\ m \ r s 

Il4? n) |l2 = T—ja / / = T-,5 / / ^ +1 (s)/ 2 (t)^ S , 

l'' i l ./o Jt Jo Jo 

where we've changed to new variables t = \m\r, s = \m\p, and used Fubini's Theorem. We 
will now estimate the above expression in two ways to show that it goes to zero when \m\ + \n\ 
increases. 

From (A. 13) we have I%(t) < ^I n (t)I' n {t) which by integration yields: 

[ l2 n(t)dt<^ l2 n(s)- (3-3) 

Next, using t < s, (3.3), and (A. 12), we get for n^O: 

, , 1 r\m\ 1 f\m\ I 

On the other hand from (A. 14) we have I^(t) < 2I n (t)I' n (t), yielding: 

f I 2 n (t)dt<I 2 n (s). (3.4) 
Jo 

So using inequalities (3.4) and (A. 12) again, we get for n^O: 

, , 1 f\m\ i p\m\ i 

Finally, if n = 0, we notice that the recurrence relations (A. 4) imply: 

tll{t) = {th{t)I {t))' -tl*(t) < (th(t)I (t))' . 
Hence we get an integral estimate: 

Il{t)tdt < sh{s)I {s) < s/ 2 (s), 



f 

Jo 
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which we use to estimate the norm above as follows: 

i r\ m \ i 

\m\ Jq \rn\ 
For the norm of S[™' n ^ we observe that, after a change of variables, we have: 



i r\ m \ ft 

\S { ^ n) \\l = r^ / ll{s)Kl +1 {t)stdsdt = \\B& 
\ m \ Jo Jo 



(m,n) |i 2 
2- 



This shows that all of the operators are indeed Hilbert-Schmidt operators. Moreover we have 
the estimates: 



\\R^ n) \\l< const ; (3.5) 
V 1 + m 2 + n 2 

and 



|4 m ' n) ||3<— = ^ llot . (3.6) 



,(m,n) || 2 ^ con St 

y/l + m 2 + ri 2 

It follows by the remarks at the beginning of the proof that Q is compact. Thus the proof 
of the theorem is complete. 

□ 

Theorem 3.2. The operator Q, defined by (2.9), is a p-th Schatten- class operator for all 
p>3. 

Proof. Notice that the p-th Schatten norm of Q can be estimated as follows: 

\\Q\\l< const (II^IK+II^I^+EII^IK- ( 3 - 7 ) 

m,n,i,j n,i 

This is because Q is (essentially) a direct sum of two by two matrices with entries made up 
of the ten integral operators we studied above, see (3.1). 

To bound H-R^'^Hp and the other norms we use the following interpolation estimate for 
the p-th Schatten norm: if a is a Hilbert-Schmidt operator and p > 2 then 

IMIp < IMI2 IHr 2 - (3-8) 
The estimate easily follows from the definition of the p-th Schatten norm. We have already 
obtained estimates on the Hilbert-Schmidt norms of R^J 1 '^ and the other operators in (3.5), 
(3.6), and (3.2), so by the above interpolation we need estimates on the operator norms. 
The main tool used to establish such estimates for the operator norms of integral operators 
is the Schur- Young inequality, see [7]. It is stated in the lemma below. 

Lemma 3.3. (Schur-Young Inequality) Let K, : L 2 (Y) — > L 2 (X) be an integral operator: 



Kf{x)= I K(x,y)f(y)dy 

Then we have: 



\\JC\\ 2 < [sup / \K(x,y)\dy) (sup / \K(x,y)\dx 

\x£X Jy J \ydYJx 
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The kernels of our integral operators are products of modified Bessel functions, and the 
difficulty here is to estimate the integrals of such products. The main technical step in those 
estimates is summarized in the following lemma. 

Lemma 3.4. Consider the following expressions for m^O: 

l[ m ^ = -L sup [ K n+1 (s)I n (t)tdt, 

|^|0<s<|m| JO 



:\m\ , 

r\ m \ 

j(m,n) = __ / K n+l (s)I n (t)sds. 



r-(m,n) 1 

\ m \ 0<t<\m\ Jt 

There is a constant such that for i = 1,2: 



j[m,n) < const 



yl + m 2 + n 2 

The proof of this lemma will be postponed until we finish the main line of the argument. 
Now we turn to estimating H-R^'"^!! for m/fl. Using Lemma 3.3 we have: 



I E)( m > n ) II 2 

II 



<m 2 (sup f K n+i (\m\p)I n+j (\m\r)pdp J ( sup f K n+i (\m\p)I n+j (\m\r)rdr 

\0<r<lJr J \0<p<lJo 



Changing variables in both integrals we get: 



sup / K n+i (s)I n+j (t)sds ) \-^— sup / K n+i (s)I n+j (t)tdt 

\ m \0<t<\m\ Jt J \\ m \0<s<\m\ Jo 

By the monotonicity (A. 9) the right hand side is biggest when % — 1 and j — 0, and so 



|^(m,n) || 2 < - 



It follows from Lemma 3.4 that: 



I D (m,n) II 2 ^ rj^m,n) sj-(m,n) 

\ K ij II S -Li ■ -L 2 



K"|| < „ „■ (3.9) 



VT + m 2 + n 2 



We turn our attention to estimating ||<S^' Tl ' r ^||. By Lemma 3.3 we have: 



\S { ™ ,n) \\ 2 < m 2 ( sup / I n+i {\m\p)K n+j (\m\r)pdp J ( sup / I n+i (\m\p)K n+j (\m\r)rdr 

\0<r<lJo J \0<p<ljp 



Clearly the expression on the right hand side of the above inequality is the same as the 
expression in the estimate of ||i?-™' n ' ) || 2 . It follows that 

1 1 || < const 
13 ~ VT + m 2 + n 2 
Consider now the case m — 0. We use Lemma 3.3 once again to compute: 

iiif ">r < (w r (*)"%) (-p /' (£)"*) < ^* 

\o<r<ij Vry / \o<p<iJ P J 1 + n 2 

and similarly 
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\0<r<l 

Either way we have for i — 1, 2: 



o<p<iio VP/ / 1 + n 2 



l( o, n ) r . const 



Combining (3.5), (3.9), and using (3.8) we get: 

Jm,n) lip - Const 



*j up — /i , 9 . 9\£^i ' 



T + m 2 + n 2 ) 2 

and exactly the same estimates for ||S'-j n ' n ' ) ||£ and Hzf 3 '™^^. Consequently, by (3.7) we get: 

m,n I 1 + m + ™ ) 2 

where the series is summable when p > 3. This concludes the proof of the theorem. □ 
Proof, (of Lemma 3.4) 

Since both K n (z) and I n (z) are symmetric for positive and negative n, see (A. 2), we will 
only need to consider the case when n > 0. 

Using (A. 14) and integrating by parts we have: 

PS PS PS 

/ I n (t)tdt < 2 / I' n+1 (t)tdt = 2sl n+1 (s) - 2 / I n+1 (t)dt < 2sl n+1 (s). 
Jo Jo Jo 

Consequently we get: 

1 f s 1 

x {m,n) = — gup / x n+l (s)I n (t)tdt < - — | sup 2sK n+1 (s)I n+1 (s). 

\ m \0<s<\m\Jo l m l 0<s<|m 

Now we bound xj m '™' 1 in two different ways. First we observe: 



Xt" l) < 2 sup K n+1 (s)I n+1 (s) < 

0<s<|m| n + L 



by (A. 11). On the other hand we have: 

x {m,n) < 2 ^ sK n+1 (s)I n+1 (s) < —-, 
|m|o< s <oo \m\ 

by inequality (A. 12). It follows that X^ m '™^ < const/Vl + m 2 + n 2 . 

We estimate X^ 71 ' 71 ^ in the same fashion, however the process is somewhat more compli- 
cated. Using (A. 16) and integrating by parts we get: 

p\rn\ p\m\ p\rn\ 

I sK n+1 (s)ds < -2 / sK' n (s)ds < 2tK n (t) + / K n (s)ds. 
Jt Jt Jt 

Using (A. 16) again yields: 

/\rn\ 
sK n+1 (s)ds < 2tK n (t) + AK n {t). 
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It follows that: 



. 1 f\m\ i 

x (m,n) = / Kn+1 ^J n ^ sds < sup (2tK n (t)I n (t) + AK n (t)I n (t)) . 



m \0<t<\m\Jt \ m \o<t<\m\ 

If n > we estimate the above expression in two ways using (A. 12) and (A. 11). First we 
have: 

x {ra,n) < 1 ( 2 | m | + 4 ) _L . 



\m\ In 



Secondly: 



If n = we have: 



x (m,n) <J_f 2+ ±\ 

\m\ \ 2n J 

|^|o<t<oo Jt 



and we need to show that the function Io(t) f°° K\(s)sds is bounded. It follows from 
the asymptotic behavior (A. 7) and (A. 8) that the limit of Io(t) J*.°° Ki(s)sds at t = is 
J °° Ki(s)sds < oo. On the other hand using L'Hospital's rule we get: 

lim J t / X 1 (s)sds = lim ov . ^ ; = hm ov ; - , w = -, 

by (A. 8) again. Thus, in a similar fashion to l[ m ' n \ we have X^" 1 '™'* < const/Vl + m 2 + n 2 . 
Therefore the proof of the lemma is complete. □ 

In conclusion we would like to remark that a somewhat more complicated proof of Lemma 
3.4 is possible without the use of the non-elementary inequality (A. 10). Estimating along the 
lines of the Hilbert-Schmidt norm bound in the proof of Theorem 3.1 it is enough to employ 
the inequalities (A. 13), (A. 14), (A. 16), and (A. 15) instead of monotonicity of K n (t)I n (t). 

Appendix A 

This section contains all of the relevant information on the modified Bessel functions. The 
references that are used are [1], [3], [11], and [13]. A short argument will be given for those 
results that are not from any of these references. 

A.l. Basic properties. The main reference of this subsection is [1]. 

The modified Bessel functions of integer order n can be defined by the following expres- 
sions: 



I n (t) = - f e tcosa cos(na) da 
n Jo 



and 



'X) 



K n (t)= I e- tcosha cosh(na) da (A.l) 



o 



where in both formulas t is a positive real number. 
Both functions are symmetric in n: 
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/„(*) = I - n {t) and K n {t) = K. n (t). (A.2) 

Consequently, without the loss of generality, it will be assumed that n is a non-negative 
integer. 

We have the following power series representation for I n (t): 

It follows that both modified Bessel functions are positive. 

They are two independent solutions of the second-order differential equation: 

d 2 x 1 dx / n 2 \ , . . . „. 

dP + m-{ 1+ ¥) x V = (A - 3) 

which is called the modified Bessel equation. 

They satisfy the recurrence relations with derivatives: 

77 77 

I' n (t) = I n+1 (t) + -I n (t) and K' n (t) = -K n+1 (t) + -K n (t), (A.4) 

as well as: 

71 ft 

I' n (t) = I n ^(t) - -/„(*) and K' n (t) = -K n ^(t) - -K n (t), (A.5) 
The Wronskian of the two functions is: 

W(K n (t),I n (t)) = det ( ^ ) = / B (t)/f n+ i(t) +/ n+1 (0^n(i) = (A.6) 

They have the following expansions near zero for n > 0: 



1 ftY , r , / n f -In (^) -7 if n = 
where 7 is the Euler-Mascheroni constant. The expansions at infinity are: 



I n (t) ~ -|= and tf B (f) ~ e"*^. (A.8) 

In the following subsections we collect less known results about the modified Bessel func- 
tions. 

A.2. Monotonicity. The modified Bessel functions have simple monotonicity properties in 
the argument t: I' n {t) > and K' n (t) < on (0, 00), which says that I n (t) is increasing and 
K n (t) is decreasing. The first inequality follows from (A.4). The second inequality follows 
immediately from the integral representation (A.l). 

Additionally there are the following monotonicity properties in the order 77: 

I n+ i(t) < I n {t) and K n {t) < K n+1 {t). (A.9) 
The first inequality was proven in [13]. It also follows from Turan - type inequality [3]: 

I n -l(t)I n+1 (t)-%(t)<0. 
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For the second inequality we estimate 

POO 

K n+ i(t) = / e ~* cosha (cosh(na) cosh a + sinh(no;) sinha) da > 
Jo 

> / e -* coshQ cosh(na)coshcKfo > / e~ tcosha cosh(na)da = K n (t). 
Jo Jo 

We also have monotonicity of the product: 

(K n (t)I n (t))' < (A.10) 
i.e. K n (t)I n (t) is a decreasing function of t, see [11]. 

A. 3. Product estimates. For n > 1 we have: 

]xmK n (t)I n (t) = 

This is a simple consequence of the asymptotics of I n {t) and K n (t) as t — > 0, see (A. 7). Since 
I n (t)K n (t) is decreasing on (0, oo), we have: 



K n (t)I n (t) < -L. (A.ll) 
In 



Additionally we have: 



tK n {t)I n {t) < tK n+1 (t)I n (t) < 1. (A.12) 

The inequality follows from (A. 9) and from the Wronskian formula (A. 6) since both terms 
on the left-hand side of that equation are positive. 

A. 4. Derivative estimates. The proofs in this paper use the following two inequalities 
with derivatives of the modified Bessel functions of the first kind. They are, for n > 0: 

In(t) < -I' n (t), (A.13) 

and 

I n - 1 <I n {t)<2I' n (t). (A.14) 

To prove them we notice that from (A. 4) we get I' n {t) — jl n (t) > 0, which gives (A.13). 
Secondly, (A. 4) and (A. 5) give: 

2I' n (t) = I n+1 (t) + I n - 1 (t) > / n _x(*) > J n (t), 

which is (A.14). 

Both inequalities above are also a direct consequence of the following stronger result of 
[11]: ^ 



I n (t) V n + 1 

For the modified Bessel functions of the second kind we have the following useful result 
from [11]: 

sK'(s) rr, k 

— < -Vs 2 + n 2 . 
KJs) 
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An analog of (A. 13) and obtainable in the same way from (A. 5) is: 

K n {s) < --K(s). (A.15) 

n 

However for our applications we only need the following estimate: combining (A. 4) and (A. 5) 
gives: 

-2K' n (t) = K n+1 (t) + K n ^{t) > K n+1 (t) > K n (t). (A.16) 
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